Interference of mesoscopic particles: quantum-classical transition 
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We analyze the double slit interference of a mesoscopic particle. We calculate the visibility of the 
interference pattern, introduce a characteristic temperature that defines the onset to decoherence 
and scrutinize the conditions that must be satisfied for an interference experiment to be possible. 
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Introduction - Interference is one of the most charac- 
teristic traits of quantum systems. As Dirac clarified 0, 
this phenomenon is rooted in the superposition principle, 
according to which different states of a single quantum 
mechanical particle interfere with each other. The sim- 
plest case is that of two states: double slit interference 
has been observed with photons, electrons, neutrons, 
atoms and small molecules 0, and recently even with 
large molecular clusters 0|. Our comprehension of the 
quantum mechanical world has been shaped, to a large 
extent, by the ideas that motivated these experiments. 

Quantum particles interfere, but classical particles do 
not, and it is not easy to understand where the borderline 
has to be placed. The size of the interfering system plays 
an important role, but it is certainly not the only rele- 
vant variable: for example, double slit interference has 
been observed with molecules, but not with protons. In 
this Letter we will analyze the interference of mesoscopic 
systems, endowed with an internal structure which leads 
to entanglement with their environment (e.g., via pho- 
ton emission). An example is a fullerene molecule flying 
between a diffraction grating and a detector. Some clas- 
sical features are apparent for such mesoscopic systems, 
yet their ability to interfere is preserved, at least to some 
extent, and can be viewed as a quantum signature. The 
main objective of this Letter will be to understand under 
which conditions a "large" system interferes and which 
ones of its dynamical variables can interfere. 

Double slit interference in a Poissonian environment - 
Let us consider a mesoscopic quantum system (molecule) , 
whose center of mass is described by a (double) wave 
packet \ipo), emerging from a double slit. The wave 
packet travels along direction +z; the slits are parallel 
to y and are separated by a distance d, along direction x 
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where p x is the x component of the momentum operator 
and IV'siit) the state emerging from one slit. We assume 



that = 0, so that |^o) is normalized. During 

its travel to the screen, the molecule emits photons and 
recoils accordingly. The internal state of the molecule 
together with the photon field plays therefore the role of 
environment: such an environment disturbs the motion 
of the center of mass, via scattering processes (typically 
photon emissions, yielding momentum kicks). We shall 
assume that the internal temperature of the mesoscopic 
system is much higher than the temperature of the pho- 
ton field. 

Let the molecule undergo momentum kicks (Api,i S 
Z) due to photon emissions. The Hamiltonian describing 
the evolution of the x component of the center of mass 
in the presence of random kicks Apk at times tk reads 
(henceforth, for simplicity, p x = p) 
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where {tk}kei, is a shot noise with density A and 
the momentum jumps Ap^ are independent identically 
distributed random variables with probability density 
W(Apk). The process is the time derivative of a 
compound Poisson process Q . Both A and W are func- 
tions of the state of the environment (for example its 
temperature of T). 

The time evolution of a wave packet which emerges at 
time to = from the slits reads 
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U^(t) being the unitary evolution engendered by the 
Hamiltonian J3J 



U 6 (t) = Texp^-ij^ dsH 6 (s)\ 
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where At k = tfc+i — tk, ApW = Y^j=i ^Pj 1S tne to- 
tal effect of k momentum jumps (Ap^ — 0), the total 
number of collisions n is a Poisson random variable with 
mean At, and T is the time-ordering operator, forcing 
earlier times (lower k) at the right. In the third equal- 
ity we used the commutation relation e~ lp At / 2m e raA P = 
e ixAp e -i{p-Ap) At/2m^ j n orc [ er mov e all kick operators 
to the far left side. 

From Q and © one gets 
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where if>o(p) = (p\ipo) and 
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is a quadratic polynomial in p with quadratic term 
p 2 Atk/2m = p 2 t/2m. It can be rewritten as 



<j)(p) = <j)(p) + — {p-p) 2 , 
2m 
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p being the value of the momentum at the extremal 
(j)'(p) = 0, that is 
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where Cfc — 1 ~ £&A characterize the emissions between 
the grating and screen. A peculiarity of this analysis is 
the presence of the same Poisson process on both branch 
waves: for an external environment, one should have con- 
sidered two independent Poisson processes, one for each 
branch wave. 

Equation JJjJ represents the convolution of the initial 
momentum wave packet with a Gaussian 
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whose spread (m/i) 1//2 becomes narrower as time t in- 
creases. For t — > oo lfTU|l reads 



ip(x,t) 
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This approximation is valid for t 3> ti|V'o(p)/V , o(p)I an( l 
implies that ipo(p) can be represented by the constant 
value V'o(p) on the screen. In the following we will always 
suppose that such condition holds (far field interference 
pattern). The interference pattern reads 



I(x,t) = (\ip(x,t)\ 
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where (• • •) denotes the average over the process 
With the initial state QJ, tpo(p) = \/2^siit(p) cos(pd/2), 
with Vsiit(p) = (plV'siit), and the far-field condition is 
satisfied for t md 2 . Under this condition, the intensity 
at the screen reads 

I(x, t) = (-J Vsiit (— ) [1 + (cos (pd))] , (13) 

where we approximated | V'siit (p) | 2 — |V'siit(wx/t)| , for 
weak enough kicks. The corresponding visibility V is 



I(x,t) = Io(x,t) 
V = IF!, a 
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F = ^exp (idJ^CkApA \ = Ve», (14) 

where Io(x,t) — (m/t)\ip s ii t {mx/t)\ 2 . In order to calcu- 
late the visibility, the features of the average (• • •) must 
be expressed in terms of the distribution of the momen- 
tum jumps W(Apk) and the Poisson times {tk}- 

By performing first the average over W(Apk) we get 
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where g(t k ) = f (d^-^-) and 

f(x) = (exp (ixAp)) Ap = / du W{u) exp(ixu), (16) 



(• • •)« denoting the average over the shot noise with den- 
sity A. This is easily computed 
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and yields the visibility (Re denotes the real part) 

V = exp(-ACt), C = / [1 - Re /(sd)] ds. 
Jo 
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By using the definition JTSJ), the "geometrical" factor- 
reads C = (1 — sinc(dAp)) Ap , where sincx = sinx/x, 
and the visibility (|18|l can be given the useful expression 



V = exp (-At (1 - sinc(dAp)) Ap 



(19) 



Notice that V < 1, because sine a; < 1. Moreover, if 
the jumps are symmetrically distributed, i.e. W(Ap) = 
W{— Ap), then f(x) is a real function and one can omit 
the real part in l|18|l. so that V = F. 

Thermodynamics - In order to calculate the visibility 
(12) from Eg. (jT^T we need to evaluate the kick rate A and 
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the probability density of the momentum jumps W(Ap) . 
Planck's blackbody formula is not valid for small atomic 
clusters and needs to be generalized on two counts. One 
is the finiteness of the number of modes (freedoms), the 
other is the reduced stimulated emission. The former 
influences the high-energy part of the spectrum, the lat- 
ter the low-energy part. In addition, there are finite-size 
effects that need to be taken into account. These are 
usually dealt with heuristically. 

If the interfering cluster can be considered (almost) 
isolated during its travel to the screen, its temperature is 
in general not in equilibrium with that of the background 
radiation field. When photon absorption from the back- 
ground radiation can be neglected (which is the case in 
which we are interested), the photon emission rate reads 
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where u> is the photon frequency, cr a b s and Cy = Nk B 
are the absorption cross section and heat capacity of the 
small particle, respectively, T its temperature and N the 
number of vibrational modes (for example, N ~ 170 for 
C 60 , iV~ 200 for C 70 ). 

The total photon rate reads 



A(T) 
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where we assume that the temperature of the molecule 
does not change appreciably due to photon emission dur- 
ing the flight. We will check the validity of this assump- 
tion later. In order to compute these quantities, we have 
to determine the u> dependence of the absorption cross 
section. We shall heuristically assume the form 



CTabs(w) = aguj 1 



(22) 



where I is a positive integer and ag a real number, and 
look for the best fit. For instance, in the case of the 
fullerenes Ceo and C70, a comparison with experiment 
yields accurate fits for a 4 = 7.04 x 10~ 66 nm 2 s 4 and 
a4 = 7.79 x 10~ 66 nm 2 s 4 , respectively By plugging 
Eq. JUJ) into (j2"0"|) one gets the series 
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and integrating term by term in Eq. H21|) one gets the 
asymptotic expansion for large N 
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Typical emission rates at T — 2500 for a time of flight 
t ~ 2 ms yields 4-5 emitted infrared photons during the 



flight in the interferometer. In such case, the temperature 
of the molecule decreases by just a few percent, which 
does not affect the emission rate, and corroborates our 
initial assumption [after Eq. I|21l) ]. 

The momentum kick on the molecule after the emission 
of a photon of frequency u) has magnitude \p\ = fko/c. By 
assuming that the emission process is isotropic, the prob- 
ability density that the molecule undergoes a momentum 
jump p reads 



W (3D \p) 
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from which the one-dimensional probability density can 
be evaluated 



(26) 



W{A P ) = J d 3 P W {3D \p)S( Px - Ap) 
c f°° duj_ 

By plugging ij2T)|) into (|TB|) we get 

f{x) = W)j ^TMsinc(^), 
and from l |18| l 

V(T,d,t)=exp(-[A(T)-G(T,d)}t) 



where 
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Note that G{T,d) - A(T), for T -> 0, so that for low 
temperatures V — ► 1. In the high temperature case, on 
the contrary, V — > 0, as expected for a classical particle. 

By plugging the series lt2"3")l into l|2"9")) and integrating 
term by term one gets 
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Visibility and quantum-classical transition - By in- 
serting the expansions Q24J1 and 1)30(1 (with I — 4, 04 = 
7.79 x 10" 66 |i|) into the visibility one gets the 
graphs in Fig.^a) for a fixed distance between the slits, 
and in Fig.^b) for a fixed time of flight. A quantum sys- 
tem, characterized by the value V = 1, tends to display 
a classical behavior, characterized by V = 0, when the 
time of flight and/or the distance between the slits are 
increased. This quantum-classical transition takes place 
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FIG. 1: (a): Visibility VvsT (0-3000 K) and t (3-20 ms) for d = 1/Ltm. (b): Visibility V vs T (0-4000 K) and d (0.01-1 /xm) 
for t — 10 ms. The dashed white level line indicates V = 1/2 and determines the "decoherence temperature" plotted in (c): 
Jdeo vs - d (0.02-1 (im) and t (1-20 ms). We plot only the physically relevant ranges of the parameters. 



at a "decoherence temperature" T^ ec determined by the 
level curve 

V = l/2 [A(T) — G(T, d)]t — In 2 = 0. (31) 

Tdec{d,t) is plotted in Fig.^c). The transition between 
the quantum and classical behavior is very sharp, both in 
Figs.nia) andQIb), and this enables us to define the de- 
coherence temperature in a clear-cut way. These graphs 
are the central results of our analysis. 

A mesoscopic system, such as a macromolecule, can 
be attributed a temperature, in the sense of Eq. (|2U[) . by 
virtue of its large number of freedoms A. In a double slit 
interference experiment, the degree of freedom associated 
with the interfering pattern [the relevant variable being 
x, see Eq. plays a special role. We now argue that, in 
general, such an "interfering" freedom is not in the same 
thermal state as the others. 

Let the experiment last for a time t (the time of flight 
of the molecule in the interferometer) and i? xch be the 
interaction Hamiltonian responsible for the coupling be- 
tween the interfering freedom and the environment. In- 
terference can be observed if 

^ H cxch dt^ <h<kTt, (32) 

where the average ( • • -} t is taken over the initial state 
of the total system yj|. The above one is a condition 
on the exchanged action and the environmental temper- 
ature. For example, if the average in Eq. (|32|l is under- 
stood in the r.m.s. sense and -ff C xch = ~£,(t) x , like in 
©, we have (J H cxch dt) = \fXi(Ap)d = Ap tot d < h 
where Aptot (T) is the total recoil due to a momentum 
random walk [2j. In such a case the first inequality in 
(|32|l is nothing but Heisenberg's inequality and this clar- 
ifies the rationale behind it. When this condition is sat- 
isfied, the macromolecule interferes. During the interfer- 
ence experiment, energy flows between the environment 
and the interfering "colder" freedom. Such a freedom, as- 
sociated with the interfering component p x , approaches 
equilibrium (at a temperature T), via momentum — and 
energy — transfer, during the momentum random walk 



process described above: eventually, visibility vanishes 
and interference is lost when the first inequality in l|32(l 
ceases to be valid. This is a rather fast process, that 
induces a classical behavior in the (relevant interfering 
variable of the) mesoscopic system. The thermalization 
process sets in afterwards, when iJ exc h<it) ~ kTt 3> h 
instead of (|32|l . and is much slower. This is, altogether, 
a remarkable picture, that adds spell to the interfering 
features of these mesoscopic systems, as well to the many 
additional problems that must be considered 10] in order 
to get a complete picture of these phenomena. 

Acknowledgements. We thank M. Arndt, S. Kuri- 
hara, I. Ohba and A. Zeilingcr for interesting remarks. 



[1] P. A. M. Dirac, Quantum Mechanics, 4th ed. (Oxford 

University Press, London, 1958). 
[2] L. Mandel and E. Wolf, Optical Coherence and Quantum 

Optics (Cambridge University Press, Cambridge, 1995); 

H. Rauch and S. A. Werner, Neutron Interj "erometry: 

Lessons in Experimental Quantum Mechanics (Oxford 

University Press, Oxford, 2000). 
[3] M. Arndt et al, Nature 401, 680 (1999); OR. Acad. Sci. 

Paris, t.2 Srie IV, 1. 
[4] D. L. Snyder, Random Point Processes, Wiley- 

Interscience, New York 1975, Ch. 3. 
[5] N. G. Van Kampen, Stochastic Processes in Physics and 

Chemistry, Elsevier, Amsterdam 1992. 
[6] K. Hansen and E. E. B. Campbell, Phys. Rev. E 58, 5477 

( 1998). 

[7] P. F. Coheur, M. Carleer and R. Colin, J. Phys. B 29, 
4987 (1996). 

[8] P. Facchi, S. Pascazio and T. Yoneda, unpublished. Care 
must be taken in order to reproduce the features of the 
absorption cross section in the visible and infrared region, 
rather than, say, in the ultraviolet one. This is because 
for the present generation of experiments [j| decoherence 
is mainly due to close-infrared photon emission. 

[9] P. Facchi, A. Mariano and S. Pascazio, Recent Research 
Development in Physics (Transworld Research Network) 
3, 1 (2002). 

[10] R.E. Grisenti, W. Schollkopf, J.P. Toennies, G.C. 
Hegerfeldt and T. Kohler, Phys. Rev. Lett. 83, 1755 



5 



(1999). 

[11] We are always neglecting the temperature of the pho- 
ton bath, see comments after Eq. Q: strictly speaking 
our environment is made up of two parts (the other free- 
doms of the molecule + the photon bath), that are not in 



equilibrium, and T is the "local" temperature of the ro- 
tovibrational and electronic modes of the macromolecule, 
responsible for photon emission. We are not considering 
the (much) slower equilibration process of the whole en- 
vironment. 



